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Abstract 

We investigate the dipole mode in a superfluid gas of Fermi atoms trapped in a harmonic po- 
tential. According to Kohn's theorem, the frequency of this collective mode is not affected by 
an interaction between the atoms and is always equal to the trap frequency. This remarkable 
property, however, does not necessarily hold in an approximate theory. We explicitly prove that 
the Hartree-Fock-Bogoliubov generalized random phase approximation (HFB-GRPA), including a 
coupling between fluctuations in the density and Cooper channels, is consistent with both Kohn's 
theorem as well as Goldstone's theorem. This proof can be immediately extended to the strong- 
coupling superfluid theory developed by Nozieres and Schmitt-Rink (NSR), where the effect of 
superfluid fluctuations is included within the Gaussian level. As a result, the NSR-GRPA formal- 
ism can be used to study collective modes in the BCS-BEC crossover region in a manner which is 
consistent with Kohn's theorem. We also include the effect of a Feshbach resonance and a conden- 
sate of the associated molecular bound states. A detailed discussion is given of the unusual nature 
of the Kohn mode eigenfunctions in a Fermi superfluid, in the presence and absence of a Feshbach 
resonance. When the molecular bosons feel a different trap frequency from the Fermi atoms, the 
dipole frequency is shown to depend on the strength of effective interaction associated with the 
Feshbach resonance. 

PACS numbers: 03.75. Ss, 03.75.Kk, 03.75.Nt 
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I. INTRODUCTION 



Among various collective excitations in an atomic gas trapped in a harmonic potential, 
the dipole mode has the unique property that its frequency is always equal to the trap 
frequency, irrespective of the interaction between the atoms. This remarkable property was 
originally proved by Kohn in the context of the cyclotron frequency of electrons in metalsQ], 
and later it was extended to an excitation spectrum of electrons in a quantum well produced 
in ALGai-a-Asylytl. Recently, Kohn's theorem has been extensively discussed in the context 
of Bo 8 e-E,„ 8 1 legion (EEC) of trapped a t0mi c g asesQ. 

Kohn's theorem is a direct consequence of the translational invariance of particle-particle 
interaction term, U int = J2i<j u ( r i — To show this, we consider an iV-body system 
described by the Hamiltonian 

N 2 i N 

1=1 zm z 1=1 

It is easy to see that the operators P a = Y^i[ m ^a^a,i — ipa,i] ( a — x -,y-, z ) satisfy the 
commutation relations [H,P a ] = £l a Pc^\- We note that the fact that [P a , U m t\ = is 
crucial to obtain this relation. If |^o) is the ground state wave-function with energy Eq, 
the three excited states generated by P a , \^?i) a = Pa\^o) (« = x,y,z), are eigenstates with 
energies E a = Q a + E , with U- mt having no effect. These excited states are referred to 
collectively as the "Kohn modes" in the theoretical literature (the "sloshing modes" in the 
recent cold atom experimental literature). 

While Kohn's theorem is exact, it is not a trivial problem as to whether or not it holds 
when the interaction is treated approximately. For example, in the mean-field approximation 
{Umt — * Umf = YaI dr'u{vi — r')n(r'), where n(r') is the particle density), the translational 
invariance of the original (exact) U int is broken, which leads to [P a , Umf] 7^ 0. Since, in 
most cases, we cannot treat many-body effects exactly, it is an important goal to obtain 
an approximation consistent with Kohn's theorem. The nature of the Kohn mode has 

— : BBBBBH. 

al.[8] and Reidl et. al.[9] proved that, in trapped Bose gases, the Hartree-Fock random 
phase approximation (HF-RPA) is consistent with this theorem at T = and at T > 0, 
respectively. 

This problem is particularly crucial for fermion superfluidity of atoms in a trap, which 
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is a topic of great current interest, because the BCS pairing approximation (Ubcs = 
-[/*|(r)*|(r)^(r)* T (r) — >■ A(r)[^{(r)^{(r) + h.c], where \l/ CT (r) is a fermion field op- 
erator and A(r) is the Cooper-pair order parameter) is almost always used in the study 
of the BCS superfluid phase. In a trap, the r-dependent order parameter A(r) destroys 
the translational invariance of the pairing interaction. Besides this, in a trapped gas 
of Fermi atoms, we need to take into account a Feshbach resonance and the associated 
molecules jiol. \v\ . The enhanced pairing interaction mediated by the molecules can be 

la Irzl [3| in atomic Fermi gases, such as 40 K and 

bl 
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used to study superfluidity 

6 Lijl^,|20,|2l|,|22|. In addition, since the strength of this effective interaction is tuna 
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external magnetic field, one can probe the BCS-BEC crossover phenomenon 
where the superfluidity continuously changes from a BCS-type to a BEC of composite 
bosons (lil 



16, 



17, U 



the crossover regime |2 




Very recently, superfluidity both in the BEC regime |27j and in 



3 11 ] has been observed in 40 K and 6 Li. In dealing with such 



systems, one must ensure that approximate theories are consistent with Kohn's theorem. 

In this paper, we give an explicit formal proof that the Hartree-Fock-Bogoliubov gen- 
eralized random phase approximation (HFB-GRPA), including a coupling between density 
fluctuations with superfluid fluctuations, is consistent with Kohn's theorem at all temper- 
atures. We prove that this is also valid for the BCS-BEC crossover. It is also valid for 
the case when one includes fluctuations around the mean-field approximation, such as first 
done by Nozieres and Schmitt-Rink (NSR)!^ at T c . This proof can be easily extended to 
deal with systems with a Feshbach resonance. Thus our results show that the NSR-GRPA 
formalism can be safely used to study response functions and collective modes in the whole 
region of BCS-BEC crossover, without fear of the breakdown of Kohn's theorem. We use 
the same formalism to explicitly verify that a zero frequency Goldstone sound mode arises 
associated with the Bose broken symmetry. 

We also show that, as expected, Kohn's theorem is no longer valid in the presence of 
Feshbach resonance if the molecular bosons feel a different trap frequency Qb from the trap 
frequency Qp which the Fermi atoms feel. In this case, the dipole mode frequency does 
depend on the strength of the effective interaction as we go through the BCS-BEC crossover 
regime. 

We give the explicit forms for the Kohn mode eigenstates in the BCS-BEC crossover 
region, without and with a Feshbach resonance. In the case of a Feshbach resonance in the 
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crossover region, we show that the Kohn mode eigenstate is not simply a rigid center of 
mass oscillation of the static condensate profile. 

This paper is organized as follows. In Sec. II, we prove that HFB-GRPA is consistent 
with Kohn's theorem at all temperatures, in the absence of the Feshbach resonance. The 
extension of this proof to the NSR-GRPA treatment of fluctuations is also discussed. In 
Sec. Ill, we generalize the proof to include the appearance of molecules associated with a 
Feshbach resonance. We also consider the case when the Fermi atoms and quasi-molecules 
feel different trap frequencies. Throughout this paper, we set % = 1 for simplicity. 



II. KOHN'S THEOREM IN THE BCS APPROXIMATION 
A. Hartree-Fock Bogoliubov approximation 



We consider a two-component Fermi gas trapped in a harmonic potential described by 
the Hamiltonian, 

H = E/ dr ^)[~^ + Krap - fW) - Uj »{(r)*{(r)^(r)* T (r), (2.1) 

where ^(r) is a fermion field operator with pseudo-spin a =j, j, and fi is the chemical 
potential. U is an s-wave pairing interaction, and V tTap (r) is an anisotropic harmonic trap 
potential given by 

Wr)= E ~™^a- (2-2) 

a=x,y,z ^ 

The proof of Kohn's theorem explained in the introduction can be easily extended to the 
second quantized Hamiltonian in Eq. (J2.1)) . When we define 

P» = Y. [ dr^l(r)[mn a f a ~tp a }^ a (r) (a = x } y } z) } (2.3) 

this operator satisfies [H, P a ] = Q a P a . Thus, the state P a \^o) has the excitation energy Q a , 
if l^o) is the ground state wave-function. 

In the Hartree-Fock-Bogoliubov (HFB) approximation, Eq. (J2.1j) reduces to 

#hfb = E/ drVl(r)[h (r) - ^n(r)]* ff (r) + J dr[A(r)*}(r)*{(r) + h.c], (2.4) 

where ho(r) = — j^- + Vt rap (r) — /i is the Hamiltonian for non-interacting Fermi atoms. Hufb 
includes two mean fields, involving the local density of Fermi atoms n(r) = J2a(^a( r )^ o-(r)) 



as well as the off-diagonal field associated with the Cooper-pair order parameter A(r) = 
— Uji&i fr)\l/|fr^|32|. In the well-known Nambu representation, -£/hfb is conveniently written 
as|33l| 

#hfb = J dr&{r) [h (r)a 3 - jn(r)a 3 + A(r)m] (r), (2.5) 
where ^ f (r) = (*t(r), *|(r)) and 

I ^t(r)\ 

¥(r) = (2.6) 

are the two- component Nambu field operators, and a a (a = 1,2,3) are the Pauli matrices. 
In Eq. (J2.5|) . we take the order parameter A(r) real and proportional to the <7i-component. 
This choice is always possible, in the absence of supercurrents or vortices. 

The HFB Hamiltonian in Eq. ()2.5|) can be diagonalized using the solution of the well- 
known Bogoliubov de-Gennes (BdG) equations, 

h* n {r) = E n V n {r), (2.7) 

where 

h = h (r)a 3 --n(r)a 3 + A(r)a 1 . (2.8) 
The energy E n of the two-component wave-function 

/ u n (r) \ 

*»(r) = ( ) (2.9) 
Vn{r) 

can be both positive and negative. Actually, a negative energy state \& n <o {—E n < 0) is 
related to a positive energy state ^/ n >o {E n > 0) by 

*n<0 = «t 2 *„> . (2.10) 
Thus, we need only solve for the positive energy solutions of the BdG equations. 



The Bogoliubov transformation 



34j is given by 



*(r) = £ [*n(r)7n T + •V a *„(r) 7 ij, (2.11) 

E n >0 



in which case Eq. (12. 4|) can be diagonalized as 



#HFB= £ E nlLln<T- (2.12) 
E n >0,cr 



Here 7 ncr is the annihilation operator of a Bogoliubov quasi-particle. When we define the 

fermion operators 7„>o = 7nt (E n > 0) and 7„<o = 7L (E n < 0), Eqs. (j2.11|) and (|2.12jl are 
rewritten, respectively, as 

*(r) = 5> n (r) 7n , (2.13) 

n 

Hn FB = J2 E nll%, (2-14) 

n 

where the summations in Eqs. (|2.13|) and (j2.14j) are taken over all the eigenstates with both 
positive and negative energies. In the following sections, we use the Bogoliubov transforma- 
tion given by Eq. (|2.13j) . 

The two mean-fields A(r) and n(r) are determined self-consistently by 

A(r) = U - f: *t (r) ffl * n (r)[l - 2/ B ], (2.15) 

E n >0 

n(r) = -J2 *t(r)(r 8 * B (r)[l - 2f n ] + 6(0). (2.16) 

E n >0 

Here /„ is the Fermi distribution function with energy E n . As usual, we need a cutoff uj c in 
the energy-summation of the gap equation (|2.15|) . The divergent term 6(0) in Eq. ()2.16|) is 
ultimately canceled out by the divergent term involved in the first term. 



B. Generalized random phase approximation (GRPA) 



The static HFB Hamiltonian in Eq. (J2.4j) neglects fluctuations in both the density and 
Cooper channels. These fluctuation effects can be included by including interactions left 
out in the HFB approximation. To describe such interactions, it is convenient to introduce 
a generalized density operator p a = ty>(r)a a i$?(r), where a = 1,2,3 represent the amplitude 
fluctuations of the order parameter, the phase fluctuations of the order parameter, and 
density fluctuations, respectively j^. Then interactions involving fluctuations in Eq. ()2.1j) 
can be written as [13, S, [3 

= ~jj dvp a (v)p a (v) (a = 1,2,3). (2.17) 

Here Uf L and ?7j L are interactions in the Cooper channel, while the interaction in the density 
channel is given by U^ L . 

In the superfluid phase, density fluctuations couple with superfluid fluctuations through 
the Josephson effect. The generalized random phase approximation (GRPA) is a RPA kind 



of approximation which includes this additional coupling[36]. When we use the GRPA to 
treat Uf L in Eq. (I2~T7I) in a linear response calculation, we find that the response in the 
generalized density is described by 

U 3 
2 



8 Pa {r,u) = ~j^ J dv'Tll p {v,v',u)5 P p{v',u) (a = 1,2, 3). 



(2.18) 



/3=i 

Here Tl^p is the zero-th order generalized density correlation function defined by 



nyr,r» = -i^ 



dte^([p a (r,t),pp(r',0)]). 



(2.19) 



As we discuss elsewhere |17||. the ordinary density correlation function is II33, while Ilju 
and II22 describe amplitude and phase fluctuations of the order parameter, respectively. 
The off-diagonal components of Eq. (j2.19|) represent coupling between these fluctuations. 
For example, n!^ expresses a phase-density coupling, originating from the Josephson effect. 
Thus, Eq. fl2.18j) describes a collective mode in terms of the amplitude (<5pi) and phase (Sp2) 
oscillations in the Cooper-channel, as well as the ordinary density oscillation (Sps). 

In a uniform gas, where we can use 5p a (r, uj) = e* q ' r 5p Q (q, uj), Eq. ()2.18|) reduces to the 
3x3 matrix equation 

( 5pi(q,w) ^ 



l-|n°(q,o;) 



0. 



(2.20) 



Here il°(q, uj) = {n^q, uj)} (a, (3 = 1,2,3) is the 3 x 3-matrix correlation function in mo- 
mentum space (The detailed expressions for n° /3 (q, cj) are given in Ref. [17].). The solution 
of the 3 x 3 matrix equation in Eq. ()2.20|) can be shown to be the same as the poles of the (3 x 3 
matrix) GRPA density correlation function, n(q, uj) = 

In a harmonic trap, n° (3 (r, r', cj) can be calculated from the analytic continuation of the 
corresponding two-particle thermal Green's functions, 



n^(r,r ,iv n 



\ Tr [°"a^( r > r '> iuJ i + Wn)<7pG(r', r, 







(2.21) 



where iui and iv n are the usual fermion and boson Matsubara frequencies, respectively. 
G(r,r f ,iui) is the 2 x 2-matrix single-particle thermal Green's function, given by 



G{r,v',iui) = 



iui - E ri 



[2.22) 



where \Pn(r) are the Nambu eigenstate spinors defined in Eq. (|2.7|) . After doing the iuj\- 
summation in Eq. (|2.21j) . we make the usual analytic continuation iv n — > uj + = to + i5. The 
resulting diagonal correlation function is given by 

E n > — E n 



',u*.-{E n ,-E n )* { 



[f n -frA *i,(r)a a y n (r) *t(i>A'(r' 



(2.23) 



where we have used the relation between positive and negative energy states in Eq. (|2.10|) . 
Similarly, the off-diagonal correlation functions n°_^ involving coupling between fluctuations 
are given by 



n? 2 (r,r» 



E 



uj\ - (E n > - E n y 



1 to+ — {E n i — E n 



[fn - f n >]\*Ur)a 2 * n (r)] M(r')a^ n/ (v' 



(2.24) 



(2.25) 



n° 3 (r,r» 



run' i 



(E n ,-E n y 



[fn - U]\9l(*)°2*n{T)] faW^A* 



(2.26) 



/L^ , ,2 



uj + — [E n i — E n 



[fn - / n >fc(r)a 3 ^(r)l M(r')a 2 y n ,(r' 



(2.27) 



n? 3 (r,r» 



E 2 \ \v \f» - M [*t(r)^i*„(r)] [*t (r>s<M*0 

nn' W + ~ W _ ^ 



(2.28) 



IlSiCry.a;) = E a \ \ , 3 Lf»-/^][^(r)^»(r)]fct(rO^('0 



(2.29) 



When we substitute Eqs. f|2.23|) - ()2.29|) into Eq. ()2.18j) . we obtain the following coupled 
equations for the collective modes 



2 ^ <4 - (s ft / - E n y 

u + (n\a 2 5p 2 (r',uj)\n'} + {E n , - E n ) (n\a 3 S p 3 (r' , uS)\ri) 



^(r)^^^)! \{E n , - E n )(n\cT 1 5p l {r\u)\n') 

(2.30) 
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6p 2 {r,uj) = -~2^ 



#1, (r)cx 2 # n (r) -lu + (n|(7i<Jpi (r', u) \ri) 



2^,u\- (E n , - E n y 
+ (E n > - E n )(n | a 2 Sp 2 (r',uj) \ n') - uj + (n\a 3 S p 3 (r' , u))\n') 



(2.31) 



Sp 3 (r,u) 



U 



E 



fn fn' 



2^ul- (E n , - E n y 
u + (n\a 2 6p 2 (r',uj)\n'} + (E n > - E n ) (n\a 3 S p 3 (r' , u)\ri) 



(2.32) 



Here we have introduced the notation 



a a 5 Pa (v')\n') = J dr'¥ n (r r )a a 5p a (v')^ n ,(v') 



(2.33) 



In HFB-GRPA formalism, the collective modes are determined by the solutions of Eqs. 
([2.30)1 - ([2.32)1 . The single particle excitations E n , particle density n(r), and order parameter 
A(r) are obtained by solving Eqs. ([2. 7)1 . (j2.15[) and ([2.16)1 self-consistently. 



C. Goldstone mode 

Since we have assumed a short-range effective pairing interaction in Eq. ([2.1)) . a cutoff u c 
is necessary in calculating correlation functions related to superfluid fluctuations, e.g., Il^- 
A similar cutoff is needed in solving the gap equation ([2.15)1 . However, we must take some 
care in how we introduce this same cutoff u c in calculating the correlation functions. For 
this purpose, a crucial role is played by involving Goldstone's theorem, which describes the 
appearance of an excitation in the ordered state associated with a broken continuous sym- 
metry. We want our HFB-GRPA formalism to satisfy this theorem and thus it is reasonable 
to introduce u c to ensure this happens. 

In the superfluid phase, the ground state is chosen from an infinitely large number of 
degenerate candidates for ground state, that are characterized by the phase of order param- 
eter A(r)e J?i (0 < cf) < 27t). The Goldstone mode is the excitation from an assumed ground 
state (where we set = 0) to the other degenerate states (0 ^ 0)[39(. As a result, the 
Goldstone mode is physically described by a phase oscillation of the order parameter with a 
zero excitation energy. The phase fluctuations of order parameter described by Sp 2 {r, uj = 0) 
in Eq. ([2.31)1 are indeed decoupled from other fluctuations at u = 0, because one can show 
n? 2 = = n° 3 = n° 2 = (see Eqs. flE23MZ2!fl ). Taking (5 Pl ,5p 2 ,5p 3 ) = (0,5p 2 ,0), we 



9 



are left with a single equation at uj = 0, 



^,(r)<7 2 *n(r)l(n|<72<J/02(r / ,0)|7i / ). 



(2.34) 



We introduce a cutoff a> c in the energy-summation (n-summation) by attaching the step 
function Q(u c — \E n \) to the Fermi distribution functions /„ and f n i in Eq. (|2.34|) as follows; 



fn —> fn = ©(^c — \E n \)f n , 
fn' — > fn' = @(^ c — \E n i\)f n i. 



[2.35) 



To prove that this prescription is consistent with the BCS gap equation (|2.15p . we next 
show that Sp2{r, 0) = A(r) is a solution of Eq. ()2.34|) when the order parameter satisfies 
Eq. From the commutation relation [h, 03] = — 2zA(r)o~2 (where h is defined in Eq. 

we obtain 

(n\a 2 A(r')\n') = % -{n\[h, a 3 ]\n') = % -{E n - E n ,)(n\o- 3 \n'). (2.36) 
Substituting <5p 2 (r, 0) = A(r) into the RHS of Eq. fl2~MD (= S(r, 0)), we find 
5(r,0) = YE[/»'-/»][*i'( r )^*»( r )]<nk3|n'> 

nn' 

= -^E/n*Ur)^i*n(r) 



U 



s„>o 



£ *t(r)a 1 *„(r)[l-2/„]. 

B„>0 



(2.37) 



In this calculation, we have used the completeness condition, J2 n ^n( r )^n( r ') = <K r ' ~~ r )- 
From the gap equation (|2.15p . Eq. ()2.37|) is found to be precisely equal to A(r). 
As a result, we have shown that at uj = 0, 

f 5 Pl (r,0)\ ( N 



5p 2 (r,0) 
\5p 3 (r,0) J 



A(r) 

V J 



(2.38) 



which describes the Goldstone mode. We conclude that Eq. (J2.35J) is consistent choice of 
the cutoff uj c with the gap equation (|2.15jh in the sense that the HFB-GRPA leads to linear 
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response functions which satisfy Goldstone's theorem. In the following sections, we use Eq. 
(|2.35J1 in evaluating the superfluid fluctuations described by 5p\ and 5p 2 in Eqs. ()2.30|) and 

(ESU). 



D. Kohn mode in Trapped Fermi superfluids 



When the frequency of a collective mode is non-zero, density fluctuations couple with 
superfluid fluctuations, so that we have to solve the coupled equations ()2.30j) - ()2.32|) . In this 
section, we prove that 



Spi(r,Q x 



5p 2 (r,tt x ) 
\ 5p3(r,tt x ) ) 



I 



\ 



—2imQ x xA(r) 
\ ~d x n(r) J 



(2.39) 



are the solutions of Eqs. ()2.30|) - (j2.32|) . with frequency u equal to the trap frequency fl x . We 
discuss the physical meaning of these solutions after giving the proof. 

From the two commutation relations [x, h] = ^(J 3 d x and [d x , h] = cr 3 mQ x x + <Ti<9 x .A(r) — 
^a 3 d x n(r), we obtain 



(E n > — E n ){n\a 3 x\n') = 2i(n\<j2A(r)x\n') H (n\d x \n'), 



m 



[2A0) 



U 



(E n , - E n )(n\d x \ri) = mn 2 x (n\a 3 x\ri) + (n|<7i(3 x A(r))|n') - -(n\a 3 (d x n(r))\n'). (2.41) 
These two equations give 

\E n , -E n ) 2 -n 2 x ](n\d x \n'} = {E n , - E n ){n\ax{d x l^{T))\n) 

+ 2imil x (n\a2A(r)x\n) 
U 



- {E n > - E n )(n\a 3 — (d x n(r))\n), 



[2A2) 



and 



m 



(E n i — E n ) 2 — Vt 2 x (n\a 3 x\n') = (ra|criA(r)x|n) 



+ 2im(E n i — E n )(n\<j2A(r)x\ri) 
~ ( n l° r 3^-(«9^(r))|r2). 



^2.43) 
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We next substitute Eq. ()2.39j) into the RHS of (|2.30p . which we denote as Si(u). At uj = Q x , 
Si can be reduced to, by using Eq. (|2.42|1 . 

^(n.) = ^E[/»-/«'][*l( r ) <7 i*»( r )]< n l 9 *'l n ') 

nn' 

= -^E/4(^( r ))^n(r)+^(r) ( r 1 5A(r)" 

* n 

z n 
r r t<; c 

= E [*Ur)<Tl*n(r)/„ + ^(^(72(7^2^(^(1 -/ n ; 



E7 



d* T E *i(r)ai*„(r)[l-2/ B 



d x A(r) 



<Jpi(r,fi x )). 



^2.44) 



Thus our trial solution in Eq. ()2.39|) satisfies Eq. ()2.3U|) for to — fl x . 

In the same way, we also find that Eq. (|2.39|) is a solution of Eq. ()2.32|) . Indeed, 
substituting Eq. (I2~l2l into the RHS of Eq. (l2~32l with lo = Q x (= S^Qa)), we obtain 

S 3 (P X ) = ^Et/n- /n'][^(r)a 3 ^„(r)](n|^|n') 

nn' 

= -^E/n*nW^n(r) 
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9 X E *n(r)^3*n(r)(l-2/ n ) 



S n >0 



-—d x n(r) (= Sp 3 (r,Q x )). 



(2.45) 



The RHS of Eq. (|2~3TJ) (= S^Ur)) can be evaluated by using Eq. (I2~l3l . We find, for 

W = Or; 



— mOsEt/" ~~ /*»'] *n'( r ) cr 2^n(r) (n|(7 3 x|n') 



X 



-iUmQ x E *n(r)<7 1 * n (r)(l-2/ n ) 

E n >0 

-2imfl x A(r)x (= 5 p 2 {r , Q x )) . 



(2.46) 



From the results in Eqs. ()2.44|) - ()2.46|) . we have shown explicitly that fluctuations as given 
in Eq. ()2.39|) are indeed a solution of the coupled equations (|2.30|) - ([2.32|) with frequency 
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Q x . This is the Kohn mode. Our calculation shows that in the BCS pairing approximation, 
the Kohn mode at the trap frequency is described as a coupled oscillation of density, phase, 
and amplitude of order parameter. We note that this contrasts with the Goldstone mode at 
uj = 0, which is associated with a pure phase oscillation. 

In the study of BEC of trapped Bose gases, it has been shown that the Kohn mode is the 
rigid center of mass oscillation of the static condensate and non-condensate distributions (| . 
We now discuss the physics of the Kohn mode in Fermi superfluids, as described by Eq. 
()2.39|) . In fermion superfluidity, the total local density profile at t is given by 

n(r, *) = n(r) + CRe[5p 3 (r, n x )e in *% (2.47) 

where C is a constant determining the amplitude of the oscillation. Substitute the third 
component of Eq. 1)2.39)1 into Eq. (J2.47)) . we obtain 

n(r, t) = n(r) — C— d x n(v) cos(Vl x t) 

~ n{r-e x C— cos(Q x t)). (2.48) 

Here e x is the unit vector in the x-direction. The Kohn mode solution in Fermi superfluids 
is thus found to be the center of mass motion of the total density, just in the case in Bose 
gasesp]. Similarly, the oscillation of the order parameter is given by 

A(r,t) = A(r) + C5A(r,t). (2.49) 

The second term is related the amplitude (dpi) and phase oscillations (5p2) of the order 
parameter, given by 

SA(r,t) = -^[Re[6 Pl (r,n x )e^ t ]- i Re[6p 2 (r,n x )e l ^ t ]]. (2.50) 

Using the first and second components in Eq. ()2.39)) . we find (working to first order in the 
fluctuations) 

A(r, t) = A(r) - C—d x A(r) cos(tt x t) + iCUmxVt x A(r) sin^t) 

~ A(r - e x C^- cos(tt x t))e iCUmxn * s ' m{nxt) . (2.51) 

Since |A(r,t)| 2 describes the density of the Cooper-pair condensate, we find that this con- 
densate oscillates in the same way as the total density profile given by Eq. ([2.48)1 . Equation 
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()2.51|) also shows that the Kohn mode is accompanied by a phase oscillation with the same 
frequency Q x . This is due to the Josephson effect, which couples the density oscillation with 
the phase oscillation of the Cooper-pair order parameter. 

It is easy to show that the other two Kohn modes with the trap frequencies fl y and fl z 
are, respectively, given by 

f 6 Pl (r,n y )\ I d y A(r) N 

-2imQ y yA{r) , (2.52) 



5p 2 (r, Q y ) 
\5p 3 (r,Q y ) J 



V 



u 



d y n(r) J 



5p 2 (r,n z ) = -2imVL z zA(v) ■ (2.53) 
\5p 3 (r,Q z ) ) \ -^d z n(r) J 

The above proof can be immediately extended to the strong-coupling theory developed 
by Nozieres and Schmitt-Rink[24] at T c in the context of superconductivity. In the NSR 
theory, fluctuations around the mean-field order parameter are taken into account within 
t he Gaussian — Q . One solves the gap equation together with the equation 
for the number of particles, which includes the effect of fluctuations around the mean field 
approximation. While the gap equation has the same form as Eq. (I2.15|) in the NSR theory, 
the chemical potential p can be very different from the Fermi energy in the strong-coupling 
regime. When we consider collective modes in NSR-GRPA theory, we again obtain the linear 
response equations in Eqs. ()2.30|) - ()2.32j) . The only difference is that the chemical potential is 
now determined by the equation for the number of particles, including the fluctuation effect 
in the Cooper channel. However, we note that the above proof based on Eqs. (|2.3Uj) - (j2.32j) 
is always valid irrespective of the value of the chemical potential. Thus, even if the chemical 
potential remarkably deviates from the Fermi energy due to the strong-coupling effect in the 
NSR-GRPA theory, we again obtain the Kohn modes with the frequencies Q a (a = x,y, z), 
as well as the zero frequency Goldstone mode. Thus NSR-GRPA linear response formalism 
is found to be consistent with both Kohn's theorem and Goldstone's theorem, and can be 
used to study the BCS-BEC crossover, without violating these important theorems. This is 
a crucial requirement of any approximate many-body calculation of response functions. 
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III. KOHN MODE IN THE PRESENCE OF A FESHBACH RESONANCE 



A. Coupled fermion-boson model 



A Feshbach resonance can 
in trapped Fermi gases 



19 



20 



)e used to produce a strong attractive (pairing) interaction 



21 



. The Feshbach resonance is associated with a dimer 



bound state, which is a boson. This boson can enhance the pairing interaction between 



atoms. To inc 
13 



14. 



15 



16 



ude this effect in a simple way, the coupled fermion-boson model 



17 



18, 



41 



10 



11 



12, 



42j is very convenient, which is given by 
V 2 _ r . _ + , „ V 2 



H = £/ «t (r)[-— + K rap - /i]*„(r) + / dr$t(r)[-— + 2u + V*. - // B ]$(r) 



c/r 



$ t (r)\I/ i (r)\I> T (r) + /i.e. 



? / dr*{(r)*J(r)^(r)* T (r) 



rti 



(3.1) 



Here $(r) is a molecular Bose field operator describing molecules associated with the Fesh- 
bach resonance. g r is the coupling constant of the Feshbach resonance, in which two Fermi 
atoms can form one Bose molecule. In turn, the molecule can dissociate into two Fermi 
atoms. 2v is the threshold energy of the Feshbach resonance; the value of 2v can be varied 
by an external magnetic field. In this model, U describes a non-resonant attractive interac- 
tion. Since one molecule consists of two Fermi atoms, we take M = 2m. We also impose the 
conservation of the total number of atoms as N = Np + 2N-B, where N-p and represent the 
number of Fermi atoms and Bose molecules, respectively. This condition has been included 
in Eq. (|3.1j) . with the chemical potentials \x and \ib = 2/i. V t f ap is a harmonic trap potential 
for molecules, given by 

(3.2) 



a=x,y,z 



If the molecules and atoms feel the same trap frequency, we have Q B = Q a . When the 
hyperfine states of atoms involved in the molecule are different from the atomic hyperfine 
states ty a (r), the molecule may feel a different trap frequency. 

Although the coupled fermion-boson model in Eq. (J3.1)) has the different form from Eqs. 
(jl.lj) and ()2.1j) due to the presence of molecular bosons, we can show that Kohn's theorem 
is exactly satisfied using Eq. ()3.1|) in the special case when Q B = Q a . The generators to 
excite the Kohn modes are now given by 



Pa = I dT^l{v)[mVt a r a -ip a ]^ a (r) 

ST J 
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+ J dT&(r)[Mn%r a -ip a ]$(T). (a = x,y,z) (3.3) 

Using Eqs. ()3.1|) and (j3.3|) . one can derive [H,P a ] = fl a P a (a = x,y,z). Thus P a \^/ ) 
(a = 1,2,3) describe excited states with frequency Q a , when |^ ) is the ground state 
wavef unction. 

The HFB Hamiltonian for the coupled fermion-boson model in Eq. (|3.1j) has the form 

HhFB = #HFB + -^HFB- (3-4) 

Here the fermion Hamiltonian H^ FB has the same form as the BCS Hamiltonian in Eq. 
()2.4j) . except that the Cooper-pair order parameter A(r) is now replaced with the composite 



order parameter, 



11 



12 



15 



17] 



A(r)=A(r)+0 r m (r), (3.5) 



consisting of the Cooper-pair A(r) plus the molecular condensate 4> m (r) = ($(r)). In the 



equilibrium state, these two order parameters are related to each other by the identity 



44] 



|A(r) = h B <p m (r), (3.6) 

where hs = — jm + + V^ p (r) — 2/x. The Bose Hamiltonian H^ FB is given by 

#!fb = / dr5&(r)h B 5<5>(r), (3.7) 

where 5$(r) = $(r) — (3>(r)) = $(r) — <p m . The identity in Eq. ()3.6|) plays a crucial role in 
the following discussion. 

The BdG equations have the same form as Eq. (|2.7|) . except that now A(r) in Eq. ()2.8|) is 
replaced with the composite order parameter A(r). Using the solution of the BdG equations, 
the Cooper-pair order parameter A(r) and the Fermi atom density n(r) are calculated from 
Eqs. (12.15)) and (|2.16|) . respectively. The equilibrium molecular condensate </> m (r) is then 
obtained from Eq. (|3.6|) . 

The presence of molecular bosons leads to an effective interaction between Fermi atoms. 
In GRPA, this interaction introduces an additional term to the linear response equation (see 



Eq. (4.2) of Ref. [17J) 



Sp(r, u) = Spu(r, u) + <5p 9r (r, u), (3.8) 
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where 



5p= 5p 2 ■ (3.9) 

In Eq. ()3.8p . 5pu{r,uj) includes the contribution of non-resonant interaction, which has 
already appeared in Eqs. (|2.3U|) - (j2.32|) . given by 

U 



(3.10) 



Here the matrix elements of Il (r, r',u) are given by Eqs. (|2.23jl - (|2.29jl . with A(r) being 
now replaced with A(r). The effect of Feshbach resonance appears in Sp gi (r, uj) as 

6p 9l (r,u) = Q lj dr'J dr"B(r,r',u;)il (v',v",u)Sp(v",u;). (3.11) 

In contrast to the non-resonance part U in Eq. (|3.10|) . the effective interaction described by 
g^B(r, r', uj) in Eq. (J3.11j) is non-local and frequency-dependent, reflecting that it is mediated 
by Bose excitations. The interaction kernel B(r, r', uS) is a 3 x 3-matrix, where only B n , i? 12 , 
B21 and B22 are finite, indicating that it only works in the Cooper-channel. These matrix 
elements can be expressed in the form 




(3.12) 



where W is the unitary matrix, 




(3.13) 



D (r,r' ,cu) is the matrix single-particle Green's function for a free Bose gas in a trap (4^. 
Using the eigenfunctions $ n ( r ) of Kb with energy 

1- 



a=x,y,z 



n a + - 



+ 2v - 2fi, 



(3.14) 



we can write D (r,r',u) as 



L>o(r,r» = £ 



^n(r)$;(rQ 
[u + i0 + )or 3 - 



(3.15) 
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B. Kohn's theorem in the presence of a Feshbach resonance 



In solving the coupled equations in Eq. ()3.8|) for a collective mode, it is useful to note 
that the contribution from the Feshbach resonance, 5p gi (r,u), does not affect the equation 
for Sp3(r,u). Thus, we easily find from Eq. ()2.45|) that 



#Pi(r,n a 



5p 2 (r,tt x ) 
\5p 3 {r,Q x ) J 



( 



d x A(r) 



\ 



(3.16) 



— 2imtt x xA(r) 
\ -^d x n(r) J 

is a solution of the equation for 5p 3 when u = Q x . 

To see if Eq. ()3.16|) also satisfies the other two equations for 5 pi and 5p 2 in Eq. (|3.1Uj) . 
we substitute this trial solution into the first (= 5pm) and second (= Spu 2 ) component of 
Eq. ()3.1()j) . These can be evaluated in the same way as in Eqs. ()2.44|) and ()2.4(i|) . and we 
obtain at to = Q x 



Spui 

\ SpU2 



d x A(r) 



-2imVt x xA(r) 



U 
9x 



d x (h B (j) m (r)) 



U 



^ —2imQ x xh B (p m (r) J 9x 



Spu(r,Q a 



(3.17) 



Here we have used the identity in Eq. (|3.6)1 . 

Next we substitute the trial solution in Eq. (J3.16j) into the first (= Sp gi i) and sec- 
ond (= 6pg r2 ) component of Eq. (|3.11|) . and set u = Q x . Noting that the factor 
/ dr"IIo(r', r", Q x )5p(r", Q x ) can be simplified by using Eq. (|3.17|) . we find 



-g T J dr'WD (r, r', Q x )WUpu(r', Sl x ). 



(3.18) 



Substituting Eqs. (ETTHJ) . (ETToT) and (l3~TTf) into Eq. (ETTHJ) . we obtain 



5p gA (r,Q x ) = g r J2 



n Sn x 



^nd x '(h B (p m {r')) - MQ x x'h B 4> m {r' 



(3.19) 



5p 9i2 (r,Q x ) = -ig x Q 3 



$n(r) 



M£ n x'h B (f> m (r') - d x ,(h B (j) m (r')) 



(3.20) 
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The operation of hs on m (r) can be conveniently calculated by expanding the molecular 
condensate wavefunction, 

m (r) = £> n $ n (r). (3.21) 

n 

Here 3> n ( r ) are the eigenfunctions of a non- interacting Bose gas of molecules in a trap [see 
Eq. (jXTSj) ]. 

Putting all these together, we have 



e - «2 



(3.22) 



$n(r) 

£n — O^ 



Mt, n x'$ n ,{r')-d x ,<$> nl {r') 



(3.23) 



The integration over r' can be carried out by expressing x' and d x > in terms of raising and 
lowering operator of a harmonic oscillator for a molecule, 

1 



it 



2M0 B 



[-d x + MO»x] 



(3.24) 



1 



2M0 B 



[d x + MQ%x\. 



(3.25) 



Since a$ n (r) = y / n^$ n _i(r) and d^$ n (r) = + 13>n+i(r) (where rail is an abbreviation 
for {n x ± 1, r^, n^.)), Eqs. ()3.22|) and (j3.23J) reduce to, respectively, 



5p gA (r,n x ) 



<M0 B „ „ $ n (r) . , ... ,. 



'M0 B 



h [52 Vn2 



la 



+i 



n+1 



n (na;>l) 



£n O a 



(3.26) 



Sp gi2 (r, O x 



M $ n (r 

nn' *»n '-■>■ 



20 B 



(fn - O x )a + (f n + tyja* 



M 



oob Q "3t [J2 Vn x + la n+1 ^ $ n (r) + ]T 



nx«n-l 



in- 



n (n x >l) 



£ — O 



^n(r) 



(3.27) 
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When the trap frequency for molecules Q% is equal to that of Fermi atoms fl x , the factors 
£n±i/(£n =t &x) disappears because of £, n ±\ = £„, ± f2f = ± Q x . In this particular case, we 
find 



5p grl (r,n x ) 



n V ^ „ (^^x) 



'Ma 



a — a' 



g T d x (f) m (r), 



(3.28) 



<5p 9r2 (r,a 



fit J! V / ^ a "-l $ n( r 



n (n x >l) 



-2img r Q x x4> m (r) . 



(3.29) 



From Eqs. (J3.17j) . (|3.28|) and (J3.29)) . we find that the upper two components of the RHS of 
Eq. ()3.8|) are described by the composite order parameter A(r) = A(r) + g r (j) m (r) as 



(3.30) 



Spui + 5p 9l i I _ / d x A(r) 
5pu2 + dp 9l 2 J \ -2imVL x xK(r) 

Thus we have proven that Eq. (J3.16|) is a solution of the collective mode equation (J3.8)) with 
frequency u = Q x (in the special case when f2f = Q x ). We can also show that 

cm \ / ^A(r) \ 

-2imVL y yh{r) , (3.31) 
V -^d y n(r) J 



5pi(r, Q y 
5p 2 (r,Q y 









8p2 






\$P3 







\ 



-2imVt z z/S.{r) 

v ■ 



(3.32) 



V ~d z n{r) ) 

are also the solutions of Eq. (|3.8j) with frequency to = Q y and Q z , respectively. These 
explicit solutions describe the Kohn mode in the coupled Fermion-Boson model. 

We now briefly discuss the physical meaning of the Kohn mode in the presence of a 
Feshbach resonance, given by Eq. (|3.16j) . The density oscillation (see Sp 3 (r, fl x ) in Eq. (j3.16|) ) 
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is found to have the same form as Eq. ()2.48j) obtained in the absence of the Feshbach 
resonance. On the other hand, the oscillation of the Cooper-pair order parameter, 

A(r, t) = A(r) - C—d x A(r) cos{fl x t) + iCUmxQ x A(r) sin(0,t) (3.33) 

cannot be written in the form given in Eq. (j2.51|) . because the composite order parameter 
appears in the RHS of Eq. ()3.33|) . This reflects the fact that the Cooper-pair oscillations are 
strongly coupled with molecular Bose excitations through the Feshbach resonance. Indeed, 
the linear response of the Bose condensate 4> m (r) induced by the oscillation of Fermi atoms 
described by Eq. (j3.16j) is given by 

(5<Z>(r,Q x )) = || dr , D° n (r,r , ,n x )[5p 1 (r',n x ) - i5p2(r',tt x )\. (3.34) 

Here Dj^r, r', Q x ) is the diagonal component of molecular Bose Green's function defined in 
Eq. (|3.14|) . Substituting Eq. (|3.16|) (which we multiply with the factor C) and Eq. (|3.14|) 
into Eq. (|3.34|) . and using the same method used in Eqs. (|3.26j) - (|3.29|) . we obtain (when 

nf = n x ) 

(5$(r,n x )) = -C^[d x - MO^l^Afr), (3.35) 

where the molecular Hamiltonian h B is defined just before Eq. (J3.7|) . In the same way, we 
also obtain 

(^(r,^)) = -C 9 ^\d x + Mtt x x\hrih(x). (3.36) 

As a result, the oscillation of the Bose condensate associated with the Kohn mode is given 
by 

m (r, t) = <f) m {r) - C^d^k^A^)) cos(Q x t) + iCm^xgJi^A^) sm{Vt x t). (3.37) 

Putting Eqs ()3.33|) and ()3.37|) together, the oscillation of the composite order parameter 
A(r, t) = A(r, t) + g T <f) m (r, t) is given by 

A(r, t) = A(r) dJ[U + l^]A(r)) cos(Q x t) + iCmQ x x[U + 1^] A(r). (3.38) 

2 v h B ' h B 

2 2 

We note that U + describes the interaction between atoms, where involves the dynam- 

h B ' h B J 

ical effect associated with the Feshbach resonance. Indeed, in a uniform Fermi gas, when 
one neglects the kinetic energy of Bose molecules, this factor reduces to U + g T /{2u — 2/i), 
which has been previously obtained as the pairing interaction associated with the Feshbach 
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resonance 



H 12, 15, 16, |l7 



Using Eqs. ()3.6|) and (j3.21J) . we can expand the equilibrium 



composite order parameter in terms of the eigenfunctions of the Bose molecules $ n (r) in a 
trap, 

A(r) = (7 r ^a„,(l + ^n)$n(r) 

= £/3 n $ n (r). (3.39) 

n 

Substituting Eq. (EOHj) into Eq. (EOSj) . we find 

A(r,t) = 5^ / a n $„(r)-^X)A»[^ + ^]9 a! $n(r)cos(^t) 



+ tCmQ x x Pn[u + 9 f] $ n (r) sin^t) 



2 

~ £/3 ft $ n (r - e,,-[f/ + 9 f\ cos(fU))e iCm ^ ( ^ +9?/f ' l)sin(n ^ (3.40) 

n 2 £ n 7 

Equation (|3.4Uj) shows that each eigenfunction $ n ( r ) rigidly oscillates around the center 
of mass with the frequency fl x . However, since the amplitude of this oscillation in each 
component (which is given by —[[/ + f-]) is different, the oscillation of the composite order 
parameter, which is given by |A(r,t)| 2 , is not described as |A(r + e x C cos(Q x t))\ 2 . This 
is quite different from the Kohn mode in a Fermi superfluid in the absence of a Feshbach 
resonance [see Eq. (J2.51j) ]. This difference arises because the effective interaction associated 
with the Feshbach resonance involves the dynamical effect of molecular Bosons, so that 

2 

the interaction between atoms depends on energy as U + j-. In the BEC limit, where 
the composite order parameter is described by the ground state of the harmonic potential 
($o( r )), only the component with n = remains in Eq. ()3.40j) . In this limiting case, Eq. 
()3.4()|) does reduce to |A(r, t)| 2 = |A(r + e x C cos(Q x t))\ 2 . Thus in this limit, we arrive at 
the usual Kohn mode solutions for a Bose condensed gas, involving an oscillation of the 
equilibrium order parameter. 

We might note that the Kohn mode solutions which appear in the linear response func- 
tions describing the Fermi fields also appear in the Bose excitation spectrum. This is to 
be expected, since in the presence of a Feshbach resonance, one knows (ill Q that 
the collective modes of the fermions associated with Cooper-pairs are strongly coupled to 
the molecular excitations. The 2 x 2-matrix renormalized Bose Green's function D in the 
HFB-GRPA is given byQ,Q 

-i 



D(u) 



E(oj)D°(u) D°(uj). (3.41) 
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Here we have used matrix notation for the dependence on r [see Eq. (J3.11|) ]. The self-energy 
correction S(u;) includes the effect of fluctuations in the Fermi atoms, given by 



w, 



(3.42) 



where the projection operator rj[A] extracts the (11), (12), (21), and (22) components from 
a 3 x 3-matrix A. The molecular Bose excitation spectrum is determined from the poles of 
Eq. (|3.41j) . given by the condition that the determinant vanishes, 



= det 



1 - ^b°(u)W^fj 



det 


! + §[ 


U-g 2 x B{u) 


TI°(uj) 


det 


1 + f W{u) 





w 



(3.43) 

In the last expression, the determinant is taken over 3 x 3-matrices in the (p±, f>2, ps)-space. 
Equation ()3.43|) always has solutions corresponding to the poles of the Fermi linear response 
functions. Thus, the renormalized Bose Green's functions in Eq. ()3.41|) also exhibit the 
Kohn mode solutions at uj = Q a (when = Q a ). 

As discussed in Sec. II, the extension to include the strong-coupling effect based on the 
NSR theory does not destroy our proof. Thus, we can safely study collective modes in the 
BCS-BEC crossover region by using the NSR-GRPA formalism even in the presence of the 
Feshbach resonance, without any breakdown of Kohn's theorem, as long as the atomic and 
molecular trap frequencies are identical (f2^ = Q a ). 

When the molecules and atoms have different trap frequencies (Q^ ^ f2 a ), the key so- 
lutions in Eqs. (j3.28J) and ()3.29|) are no longer satisfied. In this regard, we recall that, 
in the presence of a Feshbach resonance, the dominant particles continuously change from 
unpaired Fermi atoms to molecular bosons as one goes through the strong-coupling BEC 
regime (i.e., decrease the threshold energy 2i/Q Q H |l8|, ):44j ) . As a result, the average 
trap frequency which the dominant particles feel also changes from Q a to Q^, which in turn 
must affect the frequency of the "Kohn mode" in the BCS-BEC crossover. The breakdown 
of Kohn's theorem when 7^ is not due to the approximation used in HFB-GRPA (or 
NSR-GRPA), but rather is due to the changing nature of the excitation spectrum peculiar 
to a trapped Fermi gas with a Feshbach resonance. 

Fig.l shows the calculated frequency of the dipole mode at T = in the BCS-BEC 
crossover given by an approximate theory. We find that Kohn's theorem holds well (within 
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our numerical accuracy) when Qb = &f- On the other hand, the mode frequency depends 
on the threshold energy 2v in the BCS-BEC crossover when VLb 7^ £If- It continuously 
changes from the trap frequency of Fermi atoms Vtp to that of Bose molecules fig as one 
passes through the BCS-BEC crossover regime, as one expects. 



C. Goldstone's theorem in the presence of a Feshbach resonance 



In this final subsection, we briefly discuss the zero frequency Goldstone mode in the cou- 
pled Fermion-Boson model. At u = 0, the interaction kernel B in Eq. (j3.12j) is proportional 
to the unit matrix, so that, as in the BCS model discussed in Sec. II C, we need only consider 
the phase fluctuation component 8p2 in Eq. ()3. 8 jl . This has the form 

<5p 2 (r,0) = --J rfr'n° 2 (r,r',0)(5p 2 (r',0) + || dv' J dv"D° 22 (v, r', 0)n° 2 (r', r", 0). 

(3.44) 

When we take 5p2(r, 0) = A(r), we find following the discussion in Sec. II. C that the first 
term in the RHS of this equation reduces to A(r). Using the same method used to derive 
Eq. (j3.18j) . we can integrate over r" in the second term of Eq. ()3.44j) (= Spg^) to obtain 

5p gv2 = ~9rJ rfr' J D 2 ° 2 (r,r , ,O)/i B m (r'). (3.45) 
Using Eqs. (ETToD (with u = 0) and (f3~2"T|) . we find 

nn' ^ n 
n 

= g r (f) m {r). (3.46) 

Thus the RHS of Eq. (|3.44|) has been shown explicitly to equal A(r) + g T <fi m (r) = A(r). 
This proves that 5p 2 (r, 0) = A(r) describes the zero frequency Goldstone mode. In contrast 
to the BCS model in the absence of a Feshbach resonance, however, the Goldstone mode is 
now a collective phase oscillation of the composite order parameter, including contributions 
associated with both Cooper-pairs and the molecular condensate. We conclude that the 
HFB-GRPA (and NSR-GRPA) formalism is consistent with Goldstone's theorem even in 
the presence of the Feshbach resonance. We note that Goldstone mode arises even if the 
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trap frequencies felt by Fermi atoms and Bose molecules are different. We also note that 
our result also guarantees that the renormalized Bose Green's function in Eq. (j3.41j) has a 
gapless (zero frequency) excitation, a required condition for any approximate theory used 
to study Bose condensation. 



IV. SUMMARY 



In this paper, we have proved that the HFB-GRPA is a consistent formalism with Kohn's 
theorem at all temperatures. This proof is also valid for the strong-coupling superffuid 
theory developed by Nozieres and Schmitt-Rink, as used in Ref . }l7 1 . Using the NSR-GRPA 
formalism, we can safely study the linear response dynamics of the superfluid phase in the 
BCS-BEC crossover without breakdown of this general theorem on the dipole oscillation. 
The relevance of the Kohn mode in BCS superfluids was never discussed much in the context 
of superconductivity. In the case of superfluid Fermi gases trapped in a parabolic potential, 
it is important that any approximate theory used to calculate collective modes be consistent 
with Kohn's theorem. 

We also considered, for the first time, the effect of a Feshbach resonance and the associated 
formation of molecules on the Kohn mode. When the molecules feel the same trap frequency 
as that for Fermi atoms, we explicitly proved that HFB-GRPA and NSR-GRPA lead to 
Kohn's theorem. However, when the molecular trap frequency is different from the atomic 
trap frequency (which can arise when dealing with different hyperfine states), the dipole 
mode frequency depends on the strength of the effective interaction (through 2u) associated 
with the Feshbach resonance in the BCS-BEC crossover region. This was to be expected, 
of course, since the dominant excitations continuously change from Fermi atoms to Bose 
molecules as we go through the crossover regime (i.e., decrease the base molecular threshold 
2u). This result is a clear experimental signature of different trap frequencies. 

We also have given, for the first time, a detailed discussion of the various quantities which 
are oscillating in the Kohn mode in a trapped superfluid Fermi gas, as given by our explicit 
results in Eqs. (j2.39|) and (j3.16|) . In particular, we showed that with a Feshbach resonance, 
the Kohn mode involves a much more complex oscillation [see discussion after Eq. ()3.33j) 
than without a Feshbach resonance [see discussion after Eq. (J2.47j) ]. 

The BCS-BEC crossover and the Feshbach resonance are key phenomena in current stud 
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ies on superfluidity in ultra-cold Fermi gases, such as 40 K and 6 Li. In considering these 
phenomena, one should be careful to introduce approximations which do not break Kohn's 
theorem, which is an exact result of a many-body system in a harmonic trap. Our proof on 
the consistency of HFB-GRPA and NSR-GRPA with both Kohn's theorem and Goldstone's 
theorem shows that they can be used to study the collective modes of strongly- correlated 
superfluid Fermi gases. Results of such calculations will be reported elsewhere |45||. 
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FIG. 1: Calculated frequency of the dipole mode at T = in the superfluid phase of a trapped 
Fermi gas with a Feshbach resonance. We consider an isotropic harmonic trap with Q a == $7^ 
and = Qb- The Fermi energy ef = 31.5S1 ^ is for a free Fermi gas, with N = 10, 912 atoms. We 
take U = 0.001O F , g T = 0.06ft F (which gives UN/R% = 0.35e F and gr^N/R^ 2 = 0.2e F , where 



Rp = 2e-p / mVlp is the Thomas-Fermi radius), and uj c = 161. 50^. 
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